Coulomb-induced dynamic correlations in a double nanosystem 
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Abstract 

Time-dependent transport through two capacitively coupled quantum dots is studied in the 

framework of the generahzed master equation. The Coulomb interaction is included within the 

'^ ' exact diagonalization method. Each dot is connected to two leads at different times, such that 

Q ' a steady state is established in one dot before the coupling of the other dot to its leads. By 

, ^, ' appropriately tuning the bias windows on each dot we find that in the final steady state the 

transport may be suppressed or enhanced. These two cases are explained by the redistribution 

of charge on the many-body states built on both dots. We also predict and analyze the transient 

CD ' mutual charge sensing of the dots. 
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I. INTRODUCTION 

Recent on-chip measurements show how two nearby mesoscopic conductors with httle or 
no particle exchange interact via Coulomb forces. For example a quantum point contact 
(QPC) has been used as a charge detector (electrometer) near a quantum dot (QD) [l|, |2| 
and in measurements of the counting statistics of the electrons in the dot J3|. Conversely, 
the backaction of the current flowing through the QPC on the states of the QD have been 
demonstrated |j, |5[. A ratchet effect in a serial double quantum dot (DQD) driven by the 
current in the nearby QPC has been recently reported l6|. The electrons in the serial DQD 
could also be excited by photons emitted by the QPC [7|, |8| or by phonons |9[. 

Transport experiments in a parallel DQD with tunable coupling have been performed by 
McClure et al. [lO| . Both positive and negative cross current correlations have been observed 
and related to the interdot Coulomb interaction, whereas in the noninteracting case only 
negative correlations are expected [ll| . Another effect of the Coulomb correlations in parallel 



dots is the mesoscopic Coulomb drag [12|, Il3| . Unlike the macroscopic drag effect which is 



a result of quasi-equilibrium thermal fluctuations in the drive circuit [IJ], the current in an 



unbiased dot is driven by nonequilibrium time- dependent charge in the second, biased dot 
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In the theoretical descriptions of transport in parallel DQDs each dot is coupled to two 
semi-infinite leads seen as particle reservoirs with fixed chemical potentials. When the 
lead-dot coupling is weak (tunneling regime) rate or Markovian master equations are used 
P, lloUlsMlTJI. Usually the dots are considered one- level systems and the dot-dot interaction 



is reduced to one parameter. For strong lead-dot coupling a scattering theory has been for- 



mulated 12| and also Keldysh-Green methods combined with phenomenological interaction 



|14| | or with the random-phase approximation [13|. Most of the theoretical calculations were 
performed for the steady state. 

In this paper we theoretically investigate Coulomb correlation effects in capacitively cou- 
pled parallel nanosystems both in transient and steady states regime. Conventionally we shall 
call them quantum dots, but the method we use is adaptable to any sample geometry and 
any number of leads. In our setup each QD is connected to the leads at different moments 
and due to the Coulomb interaction they mutually respond to each other's transient charging 
or discharging. The aim of this work is to describe and understand these effects. Depending 
on the initial conditions (occupations, bias voltages) the current cross-correlations may be 
positive or negative. The calculations are performed within the generalized master equation 
(GME) method for the reduced density operator (RDO) of the double dot. The formalism 



was adapted for open mesoscopic systems by several authors [18|. We used it recently to 



study the transient behavior of open noninteracting |19l-l2l| and interacting nanosystems 
J22| . The interaction is treated with the exact diagonalization method, both intradot and 
interdot, on equal footing. 



II. THEORY 

The Hamiltonian of the total system, shown in Fig. [H is 

H{t) = Hs + H{i} + Hrr{t), (1) 



where S stands for the "sample", in this case the DQD, i.e. QD^^ + QD^, and {/} 
{La, Ra, Lb, Rb} is the set of leads. i^T incorporates the sample-leads tunneling. 



n I -^ 



(2) 



with Xi{i) time-dependent functions describing the contact with the lead /. Cgi/dn are the 



creation/annihilation operators in the leads and sample respectively, and T' are model 



specific coupling coefficients.. The RDO p{t), or the "effective" statistical operator of the 
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FIG. 1: The system: A double dot and four leads. 

open sample, is defined by averaging the statistical operator of the total system over the 
states of all leads. In the lowest (quadratic) order in Ht it satisfies the GME. 
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where Ut = e~'^^^^^^^^'-^''^ is the evolution operator of the disconnected system, and p{;} 
is the statistical operator of the leads in equilibrium which is the product of the Fermi 
distributions of each lead / with chemical potential oj. Before the leads are coupled p{t) 



describes an equilibrium state of the isolated sample |19|-|22 



Hs also includes the Coulomb interaction. The interacting many-electron states (MES) 
of the isolated sample, solutions of Hs\oi) = Ea\a), are found by exact diagonalization. Each 
MES is expanded in the Fock space built on a finite number of single-electron states (SES), 
^SES- The number of electrons A^ in the sample may vary between zero and Nses and 
hence the number of MES is Nmes = 2^^'^^. Since the dots are not in direct tunneling 
contact the number of electrons in each dot, Na and Nb, respectively, are "good quantum 
numbers" for the MESs. The ground-state energies of the isolated DQD can be labeled as 
Sg{Na, Nfj). The chemical potential of a MES with N = Na + Ni, electrons, p{Na, Nb), is the 
energy cost to add one more electron to the ground-state with A^ — 1, and has to fit with 
the leads' chemical potentials in order to allow transfer of electrons. 

We solve Eq. Q numerically in the MES basis {|a)}. Using the RDO we can calculate 
the mean number of electrons and hence the charge in each dot, Qi{t), [i = a,b), and by 
taking the time derivative we obtain the currents in each lead. 



Qi{t) = '^^i^Po^n.ar.^ = Jhi{t) - JRi{t), 



(4) 



Ui 



where «„- are the MESs of the double system with Ui electrons in QDj. We can thus 
describe the partial charge and currents associated with any partition of electrons. The 
currents corresponding to each lead are identified from the last term of Eq. ([3]) 19|, |22] . A 
current is positive when fiowing from left to right and negative otherwise. 



III. RESULTS 

We use a lattice model for our system, each QD being a chain of four sites. The electrons 
are distributed on a 2 x 4 lattice, but the hopping between the chains is forbidden. The 
coupling coefficients are 

Tin = Vor,i{0)Uk) (5) 

where ipgi and 0„ being the single-particle wave functions in the leads and in the sample, 
evaluated at the contact sites labeled as and ii, respectively [19|. The parameter Vq gives 
the coupling strength. All electrons on this lattice interact with pairwise Coulomb potentials 
U/djk with djk the distance between electrons j and k and U a strength parameter. Coulomb 
forces are neglected in the leads. We use all 8 SES of the lattice to calculate all 256 MES, 



and the first 40 MES are sufficient to obtain convergent results from Eq. ([3]) |22|. Our 
energy unit is the hopping energy in the dots t/j, the time unit is h/ti), and the currents 
are calculated in units of et£)/h. We use Vq = 1.5 and U = 1. The ground-state energies 
for the DQD are ^^(0,0) = 0, ^^(1,0) = 2.38, ^^(l, 1) = 5.47, ^^(2,0) = 6.36, ^9(2,1) = 
10.08, Sg{3,0) = 12.37, etc., and thus yu(l, 1) = 3.09, yu(2,0) = 3.97, yu(2, 1) = 4.6. The 
dots being identical £g{Na, N^) = Sg{Ni,, Na). 

In the following cases both dots are initially empty and yULa = A^Lb, /U-Ra = yURb- QDa opens 
at ta = and after a charging period it evolves towards a steady state. In Fig. [2]^a,b) we 
show the current in QDa for two choices of the chemical potentials of the leads. In the first 
case fi{l, 1) < /iRa < /x(2, 0) < /iLa < At(2, 1), mea ning that in the steady state of QDa the 
main contributor to the current is the MES (2,0) [22]. QDb is coupled at tf, = 120 when a 
new transient period begins for both dots, after which all currents end up at equal values, 
considerably smaller than before tf,. So one can say the two dots are negatively correlated: 
The activation of one inhibits the other until they block each other. Fig. [21^a). In the second 
case. Fig. [2](b), we have instead yu(2, 0) < /iRa < /i(2, 1) < /iLa and only a very small current 
passes through QDa in the first steady state due to the Coulomb blockade. But the coupling 
of QDb now activates QDa, so the dots become positively correlated [lO|. 

To explain what is going on we show in Fig. [2](c) the population of the relevant states 
for the first case, calculated with Eq. (jl]). As long as QDb is closed one- and two-particle 
states of QDa are charging yielding a total charge up to Qa ~ 1.5, as further shown in Fig. 
M^d). Once QDb opens the electrons tunneling into it repel some charge from QDa and new 
MESs are being created like (1,1) and (2,1). Since Sg{2,0) > Sg{l, 1) the new two-particle 
ground state is (1,1) and hence the transition (2, 0) — )■ (1, 1) occurs, but also (2, 0) — )■ (2, 1). 
The later is possible because £^^(2,1) — £^^(2,0) = 3.72 is slightly below the bias window. 
Consequently the states (2,0) depopulate fast whereas the populations of the states (1,1) and 
(2, 1) (and (1, 2) as well) increase, as seen in Fig. [2](c). In the steady state the bias window 
is nearly empty of any MES chemical potential and consequently the currents nearly vanish. 
This is an interdot Coulomb blocking effect [10|]. The total charge in the dots converges to 
2.2 electrons. Of that 1.5 reside on two-electron states: 1.2 on the ground state (1,1), i.e. 
below the bias window, and 0.3 on excited states (1,1) and MESs (2,0). Also, about 0.6 
electrons are on three-particle states (2,1), i.e. above the bias window. 

In Fig. ^a) we show the partial currents in the leads connected to QDa carried by the 
two- and three-particle states. The former drop fast after tb during the depletion of the MES 
(2,0). Because /i(2,0) is almost in the center of the bias window jRa,2 and JLa,2 are very 
similar. The three-particle currents in QDa are more interesting. They correspond to the 
states (2,1) and, surprisingly, JLa,3 < and Jija.s > 0, meaning that QDa ejects charge in 
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FIG. 2: (Color online) (a)-(b): The total currents in the leads for two bias windows: (a) /XLa = 
/iLb = 4.25, /XRa = /URb = 3.75; (b) /iLa = /"Lb = 4.75, ^Ra = ^b = 4.35. The insets show the bias 
windows and the MES chemical potentials, (c) Partial charge for the states involved in (a) (see 
text); (d) Total charge on each dot for (a) and (b). 



both leads L^ and Ra- The resulting "lobe" shape is also seen in Fig. [2|a). The net charging 
of the (2,1) states is actually done through the leads connected to QDb, as can be seen in 
Fig. [31(b), where JLb,3 > and JRb,3 < 0, i.e. both currents flow into the dot. 

We return now to the positive correlation case. For t < i^ no chemical potential of type 
^{Na^ 0) is inside the bias window and hence no current flows in the steady state of QDa. 
The charging goes up to Qa ~ 1-9 as seen in Fig. EJ^d), with the ground state (2,0) occupied. 
When QDb is open the new states (2,1) are created and since the corresponding /i(2, 1) is 
inside the bias window they are available for transport. After the transient phase, when 
QDb is charging and QDa is discharging, all currents reach the same steady value, driven 
by the states (2,1) and (1,2) which end up equally populated. The currents in the steady 
state have now two- and three-particle components. These partial currents, shown in Fig. 
Wis), have another curious behavior. Both JLa,3 and JRa,3 are positive in the steady state, 
whereas JLa,2 and JRa,2 are negative, the net result being the total, positive current. This 
means that three-particle currents flow from left to right, but the two-particle currents go 
from right to left. The reason is that in our model the electrons are created or annihilated 
one at a time. A MES (2,1) is formed by creating one more electron to the ground MES 
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FIG. 3: (Color online) The partial currents carried by the two- and three-particle states in QDa 
and QDb: (a-b) corresponding to Fig. ^a); (c-d) corresponding to Fig. ^h). 



(1,1), and so the positive (2,1) and (1,2) currents deplete the (1,1) states. But yu(l, 1) is 
below the bias window and so the (1,1) states have to be backfed by a negative, two-particle 
current. The single particle states are not occupied and do not contribute to transport. The 
current in the circuit a is carried by MES (2,1), but not by (1,2), and the other way round in 
the circuit b. The electrons tunneling from the left leads a and b thus compete each other to 
access the (1,1) MES. In turn, when electrons leave a dot the remaining two-particle MES 
has the lowest energy Sg{l, 1) and not £^^(2,0) which is higher. Such transitions are called 



"U-sensitive processes" in Ref. llO 



The partial currents of states (2,1) and (1,2) are shown in Fig. Ill^a,b). In this case we 
use the same chemical potentials as in Fig. [2]^b), but now QDb contains one electron in the 
ground state at t = 0. After t = th QDb absorbs more charge and the double system evolves 
toward the same steady state as before. Fig. ID^c). But prior to U, although isolated, the 
initial electron is being excited by the charging of QDa. This can be seen in the Fig. Hl^d) 
where the populations of the ground state (2,1) and of the MESs containing the excited 
state of the electron in QDb denoted as (2, 1^) are displayed. The currents in the circuit a 
feel the initial electron in QDb, but also the excited states of it. Indeed the MESs (2, l^.) 
decay while the system approaches the steady state. 

Next we keep /iLb fixed and decrease /iRb relatively to the setup of Fig. M^a) , for increasing 
the bias eVb = [lu, ~ /^Rb- Fig. [5t^a) shows that the splitting of the two currents in QDa 
during the second transient phase decreases. The final current increases with VJ,, but it is 
still smaller than the steady value before th- The effect of increasing Vb on the final currents 
occurs in two steps. First the states (2,0) and (0,2) become slightly populated and tunneling 
of one electron creates three-particle currents. Fig. ^a) with /iRb = 3, 25. Then, the bias 
window approaches /i(l, 1) and eventually includes it, and tunneling on MES (1,1) amplifies 
the three-particle currents. Since VJ, is acting directly on QDb the final currents in the b 
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FIG. 4: (Color online) (a-b) The contributions of the ground state (2,1) and of the excited states 
(2, 1^.) to the currents in the leads. //La = ^Lb = 4.75, //Ra = MRb = 4.35. QDb initially contains 
one electron, (c) The evolution of total charge in each dot. (d) Populations of the three-particle 
states. 
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FIG. 5: (Color online) The effect of the bias applied on QDb on the currents JLa.Ra- (a) /^La = 
//Lb = 4.25, /iRa = 3.75. (b) /XLa = /^Lb = 4.75, /iRa = 4.35. The same line type is used for the left 
and right currents. 



circuit are larger than in a (not shown). Fig. [5t^b) shows the result of increasing V^ starting 
with the setup of Fig. EJ^b). The Coulomb blockade on QDa is still lifted when the bias on 
QDb increases. The discharging in the a arm may be large enough to produce a negative 
left current. 

Finally, one comment on the interdot distance. The interdot Coulomb interaction de- 
creases with the distance between the dots, but for simplicity we kept it equal to the lattice 
constant. Increasing it the MESs change, but similar effects were obtained by appropriately 
tuning the chemical potentials of the leads. 



IV. CONCLUSIONS 

In conclusion, we discussed time dependent charge sensing effects and computed mutually 
sensitive currents in parallel quantum dots. A steady-state transport regime of one dot 
is suppressed after connecting the second one. Conversely, the current through one dot 
increases if the charging of the second dot opens new many-body channels within the bias 
window. In particular, we predict that the transient current in the leads attached to the first 
dot may change sign when the second dot is connected. This effect can be experimentally 
tested. 

The RDO of the coupled system and the GME describe its entangled dynamics by treating 
all electrons equally. The Coulomb effects are fully included and the charging and discharg- 
ing energies are present in the MES structure. The classical charging and the quantum 
correlations are treated together. The exact many-body states fit naturally with the Fock 
space formulation of the GME. The access to individual MES allows a better understanding 
of the Coulomb-induced effects on the total currents. 
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